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Hybrid Control of a Bioreactor with Quantized 
Measurements: Extended Version 

Francis Mairet and Jean-Luc Gouze 


Abstract —We consider the problem of global stabilization of an 
unstable bioreactor model (e.g. for anaerobic digestion), when the 
measurements are discrete and in finite number (“quantized”), 
with control of the dilution rate. The model is a differential 
system with two variables, and the output is the biomass growth. 
The measurements define regions in the state space, and they 
can be perfect or uncertain (i.e. without or with overlaps). 
We show that, under appropriate assumptions, a quantized 
control may lead to global stabilization: trajectories have to 
follow some transitions between the regions, until the final 
region where they converge toward the reference equilibrium. 
On the boundary between regions, the solutions are defined as a 
Filippov differential inclusion. If the assumptions are not fulfilled, 
sliding modes may appear, and the transition graphs are not 
deterministic. 

Index Terms —Hybrid systems, bioreactor, differential inclu¬ 
sions, quantized output, process control 


I. Introduction 

Classical control methods are often based on the complete 
knowledge of some outputs y(t) of the system Il22l . By com¬ 
plete, we mean that any output yi is a real number, possibly 
measured with some noise 8i. The control is then built with 
this (noisy) measurement. These tools have been successfully 
applied in many domains of science and engineering, e.g. in 
the domains of biosystems and bioreactors ED- However, in 
these domains, detailed quantitative measurements are often 
difficult, too expensive, or even impossible. A striking example 
is the measurements of gene expression by DNA-chips, giving 
only a Boolean measure equal to on (gene is expressed) 
or off (not expressed). In the domain of bioprocesses, it 
frequently happens that only a limited number or level of 
measurements are available (e.g. low, high, very high ...) 
because the devices only give a discretized semi-quantitative 
or qualitative measurement The measure may also be 
quantized by some physical device (as the time given by an 
analogical clock), and give as a result some number among a 
finite collection. 

For this case of quantized outputs, the problem of control 
has to be considered in a non-classical way: the control 
cannot be a function of the full continuous state variables 
anymore, and most likely will change only when the quantized 
measurement changes. Moreover, the control itself could be 
quantized, due to physical device limitations. 
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The above framework has been considered by numerous 
works, having their own specificity: quantized output and 
control with adjustable “zoom” and different time protocols, 
ll20l , hybrid systems abstracting continuous ones (cf. IT71 for 
many examples of theories and applications). 

In this paper, we consider a classical problem in the field 
of bioprocesses: the stabilization of an unstable bioreactor 
model, representing for example anaerobic digestion, towards 
a working set point. Anaerobic digestion is one of the most 
employed process for (liquid) waste treatment ED- Consid¬ 
ering a simplified model with two state variables (substrate 
and biomass), the system has two stable equilibria (and an 
unstable one, with a separatrix between the two basins of 
attractions of the respective stable equilibria), one being the 
(undesirable) washout of the culture ( fl5l ). The goal is to 
globally stabilize the process toward the other locally stable 
reference equilibrium. The (classical) output is the biomass 
growth (through gaseous production), the control is the dilu¬ 
tion rate (see |23l for a review of control strategies). There 
exists many approaches based on well-accepted models (|7J), 
for scalar continuous output ED. 0 or even for delayed 
and piecewise constant measurements ED- Some original 
approaches make use of a supplementary competitor biomass 

EQ. 

In this paper, we suppose that the outputs are discrete or 
quantized: there are available in the form of finite discrete mea¬ 
surements. The precise measurements models are described 
later: roughly, the simplest one is “perfect”, without noise, 
meaning that the true value is supposed to be in one of the 
discrete measurements, and that the transitions between two 
contiguous discrete measures are perfectly known. The next 
model is an uncertain model where the discrete measurements 
may overlap, and the true value is at the intersection between 
two quantized outputs. Remark that the model of uncertainty 
is different from the interval observers approaches (|2], |fT3l ) 
for the estimation or regulation m, where some outputs or 
kinetics are not well known, but upper or lower bounds are 
known. Moreover, in the interval observer case, the variables 
are classical continuous variables. 

For this problem, the general approaches described above 
do not apply, and we have to turn to more tailored methods, 
often coming from the theory of hybrid systems, or quan¬ 
tized feedbacks (see above) ... We here develop our adapted 
“hybrid” approach. It has also some relations with the fuzzy 
modeling and control approach: see e.g. in a similar bioreactor 
process the paper IfTTI . We provide here a more analytic 
approach, and prove our results of stability with techniques 
coming from differential inclusions and hybrid systems theory 
m. Our work has some relations with theoretical qualitative 
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control techniques used for piecewise linear systems in the 
field of genetic regulatory networks (0). The approach is 
also similar to the domain approaches used in hybrid systems 
theory, where there are some (controlled) transitions between 
regions, forming a transition graph 0 , ra. 

The paper is organized as follows: Section [IT] describes the 
bioreactor model and the measurements models. The next 
section is devoted to the model analysis in open-loop (with 
a constant dilution rate). In Section |IV] we propose a control 
law and show its global stability through the analysis of the 
transitions between regions (with the help of the Filippov 
definition of solutions of differential inclusion). Finally, in 
Section E we explain how to choose the dilution rates by a 
graphical approach, and we end by giving some simulations, 
with or without uncertainty. 


we only know that y(£) is in a given range, or equivalently 
that £ is in a given region (parameter a is a positive yield 
coefficient): 

Yi ={£ eR 2 + :y.< z/(£) < yj, i=l,...,n-l, 

Y n ={£ e R+ : < y(£)}. 

where 0 = y x < y 2 < ... < y n and < y 2 < ... < y n _ v 
We will consider two cases: 

• (Al). Perfect quantized measurements: 

Vi = V i+ v V* e 1}. 

This corresponds to the case where there is no overlap 
between regions. The boundaries are perfectly defined 
and measured. 


II. Framework 

A. Model presentation 

In a perfectly mixed continuous reactor, the growth of 
biomass x limited by a substrate s can be described by the 
following system (see El, 03): 

/ s = u(t)(s in - s) - kp(s)x 
\ X = (fj,(s) — u(t))x 

where Si n is the input substrate concentration, u(t) the dilution 
rate, k the pseudo yield coefficient, and p(s) the specific 
growth rate. 

Given £ = (s, x), let us rewrite System 0 as 

£ = /(£, u(i)), where the dilution rate u(t) is the manipulated 
input. 


The specific growth rate p(s) is assumed to be a Haldane 
function (i.e. with substrate inhibition) 0: 


p(s) = p 


ks 


7 kj 


( 2 ) 


where /j, ks, ki are positive parameters. This function admits 
a maximum for a substrate concentration s = y/k$kj := s, 
and we will assume s < Si n . 


Lemma 1. The solutions of System 0 with initial conditions 
in the positive orthant are positive and bounded. 


Proof: It is easy to check that the solutions stay positive. 
Now consider z = s + kx whose derivative writes 


z = u(t)(s in - z). 

It follows that z is upper bounded by max(c(0), Si n ), and so 
is kx. Finally, if s(t) > Si n , then s(t) < 0, therefore s is 
upper bounded by max(s(0), Sj n ). ■ 

In the following, we will assume initial conditions within 
the interior of the positive orthant. 


B. Quantized measurements 

We consider that a proxy of biomass growth y(£) = ap(s)x 
is monitored (e.g. through gas production), but in a quantized 
way, in the form of a more or less qualitative measure: it can 
be levels (high, medium, low...) or discrete measures. Finally, 



o y 4 = Vi-i Vi = y i+1 y 


• (A2). Uncertain quantized measurements: 

J/. < Vi-x <y i+v Vie {2,...,n—1}. 

In this case, we have overlaps between the regions. In 
these overlaps, the measure is not deterministic, and 
may be any of the two values. 


Yi 

Y i -1 t _a _ v Yi +1 


o my ,-1 h+iVi y 

For both cases, we define (open) regular domains: 

Y z -YiMY^UYi+t), 
and (closed) switching domains: 

^i|i+i “ Yi n Yi+i 

where the measurement is undetermined, i.e. if £ £ Y\i+ \ • 
then either £ £ Yi or £ £ Y l+ i. 

For perfect measurements (Al), we have Y, = intYi, 
and the switching domains Y^ i+1 correspond to the lines 
y(£) = y i = j/ . For uncertain measurements (A2), the 
switching domains Y^ i+1 become the regions {£ £ : 

V i+1 < 2/(0 < Vi}- 

Unless otherwise specified, we consider in the following 
uncertain measurements. 


C. Quantized control 

Given the risk of washout, our objective is to design a 
feedback controller that globally stabilizes System 0 towards 
a set-point. Given that measurements are quantized, the con¬ 
troller should be defined with respect to each region: 

£(t) e Yi u(t) = Di, i = l,...,n. 


(3) 
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Here Di is the positive dilution rate in region i. This control 
scheme leads to discontinuities in the vector fields. Moreover, 
in the switching domains, the control is undetermined. Thus, 
solutions of System (0 under Control law ([3]) are defined in the 
sense of Filippov, as the solutions of the differential inclusion 

E2: 

i e H(0 

where II (f) is defined on regular domains V, as the ordinary 
function iT(£) = /(£, Di), and on switching domains Y i | i+1 
as the closed convex hull of the two vector fields in the two 
domains i and i + 1: 

H(t)=m{m,D i ),f{s,D i+ 1 )}. 

Following E). a solution of System 0 under Control 
law 0 on [0, T] is an absolutely continuous (w.r.t. 1) function 
£(t, £o) suc h that £(0, £o) = £o and £ C iT(£) for almost all 

te[0,T\. 

III. Model analysis with a constant dilution 


The asymptotic behavior of the system can be summarized as 
follows: 

Proposition 1. Consider System 0 with a constant dilution 
rate u(t) = D and initial conditions in the interior of the 
positive orthant. 

(i) If D < p(si n ), the system admits a globally exponentially 
stable equilibrium S; a (D) = (s a (D), Stn ~^ a ^ D ^ m 

(ii) If p(Si n ) < D < /i(s), the system admits two locally 
exponentially stable equilibria, a working point £ a (Z?) = 
(s a {D), s » n ~^ a (D) ^ anc i the washout = (sj n ,0), 

and a saddle point ^(D) = (st,(D), ^ see 

Figure [7] 

(iii) If D > p(s), the washout £o = (si„,0) is globally 
exponentially stable. 

Proof: See 031 . ■ 

For D £ (s)), let us define: 

qD 

y j( D ) = — [Sin — Sj(D)\, j = a, b. 



Fig. 1. Phase portrait of System m with a constant dilution rate u(t) = D£ 
n(s)) (case ii of Proposition |T}. Magenta lines: trajectories, cyan 
dashed lines: nullcline y(£) = 0, green dash dotted lines: isolines y(£) = 
y a (D) and y(£) = yi(D), red dotted line: separatrix, black vertical line: 
s = s, dark circles: stable equilibria, open circle: unstable equilibrium. 


We first consider System (Q} when a constant dilution rate D 
is applied in the whole space (i.e. u(t) = D, Vt > 0). In this 
case, the system is a classical ordinary differential equation 
which can present bistability, with a risk of washout. Let us 
denote s a {D ) and Sb(D) the two solutions, for D £ (0 ,p(s)), 
of the equation /r(s) = D, with 0 < s a (D) < s < Sb{D). For 
the Haldane growth rate defined by (0, we have: 



y a {D) and yb(D) are the growth proxy obtained respectively 
at the equilibria £ a (-D) and Cb(D) (if it exists £[ 

In order to design our control law, we need to provide some 
further properties of the system dynamics. In particular, we 
need to characterize j/(£), the time derivative of y(£) along a 
trajectory of System 0} with a constant dilution rate u(t ) = I): 

y(0 = a [D(sin - s) - kp(s)x] p!{s)x + - D)x. 


Let consider the following functions: 

p(s) - D D(s in - s) 


9d ■ si 


h J D : s i 


kp'(s) kp(s) 

D(Sin - Sj(D)) 


kp,(s) 


j = a, b, 


defined respectively on (0, s ) U (s, +oo) and (0, +oo). 

In the (s,x) plane, go (s), h a D {s) and h b D (s) represent 
respectively the nullcline y(£) = 0 and the isolines y(£) = 
y a {D) and j/(£) = yb(D) (i.e. passing through the equilibria 
£, a (D) and £j,(7J)), see Figure [I] Knowing that the nullcline 
2/(£) = 0 is tangent to the isoline y(£) = y a (D) (resp. 
t/(£) = yb(D )) at the equilibrium point £ a (£>) (resp. ^ b (D)), 
we will determine in the next lemma the relative positions of 
these curves, see Fig. Q] 


Lemma 2. Consider System (|T} with a constant dilution rate 
u(t) = D. 

(i) For s £ (0, s), we have go(s) > hf)(s): the nullcline 
y{£) = 0 is above the isoline y(£) = y a (D). 

(ii) For s £ (s, Sm )> we have gD{s) < h b D (s): the nullcline 
y(£) = 0 is below the isoline y (£) = yb{D). 


Proof: See Appendix. ■ 

This allows us to determine the monotonicity of y(f ) in a 
region of interest (for the design of the control law). 


1 if D < //.(.s, r; ,), £b(D) does not exist and yb(D) < 0. 
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Lemma 3. Consider System (JTJ with a constant dilution rate 
u(t ) = D. For £ £ such that yb(D) < y(£) < y a (D), we 
have y (£) > 0. 

Proof: See Appendix. ■ 

IV. Control with quantized measurements 

A. Control design 

Our goal is to globally stabilize the system towards an 
equilibrium with a high productivity (the productivity is the 
output afi(s)x), corresponding to a high dilution rate !),, 
where there is bistability in open loop (case (ii) of Proposition 


0}. Number n is the number of measurements (see section 
HEU and the control is such that D\ < Di < ... < D n . We 
consider the following control law, based on the quantized 
measurements y(£), and constant within a given region: 

Vf>0, Z(t)eYi 4=> u(t) = Di, (4) 

given that the following conditions are fulfilled: 

yb(Di)<y. i = l,...,n, (5) 

y a {Di)>y z i = 1,... ,n — 1, (6) 

y a (D n )>y n _ i- (7) 


These conditions make the equilibrium £, a (D n ) globally 
stable, as we will see below. In Section IV-AI we will precise 
how to choose the Dj such that these conditions hold. 

In order to prove the asymptotic behavior of System 0 under 
Control law 00. the study will be divided into three steps: 

• the dynamics in one region, 

• the transition between two regions, 

• the global dynamics. 

This approach is similar to those deducing the global dynamics 
from a “transition graph” of possible transitions between 
regions 0. 

B. Dynamics in one region with a given dilution: exit of 
domain 

We first focus on a region Y, . i < n. A constant dilution 
Di - such that Conditions (06]) for i hold - is applied. These 
conditions guarantee that the stable operating equilibrium for 
this dilution (see Proposition 0 is located in an upper region 
Yj, j > i, while the saddle point is located in a lower region 
Yfe, k < i. This allows us to establish the following lemma: 

Lemma 4. For any i £ {1...., n 1}, consider System 0 
under a constant control u(t) = Di, such that Conditions 0 
[6]) for i hold. All solutions with initial conditions in Yi leaves 
this domain, crossing the boundary y(£) = y\. 

Proof: Let us consider the function V(£) = y a (D n ) — 
j/(£) on Yj. Given Conditions 00, we get £ V,: 

y b {Di) < y. < y(0 <Vi< y a {Di). 

Since a constant dilution rate Dj is applied, we can apply 
Lemma 0 to conclude that y(£) > 0. 

Thus, V(£) is decreasing on Y, . Recalling that the trajec¬ 
tories are also bounded (Lemma 0, we can apply LaSalle 


invariance theorem m on the domain := {£ £ Yj | 
x < max(^(0), s in )/k, s < max(s(0), s in )}. Given that 
the set of all the points in Tli where V(£) = 0 is empty, 
any trajectory starting in fli will leave this region. The 
boundaries x = max(^(0), Si n )/k and s = max(s(0), Sj n ) 
are repulsive (see Proof of Lemma 0. Finally, the boundary 
y(0 = y.j corresponds to the maximum of V(£) on fli, so 
every trajectory will leave this domain, crossing the boundary 

y(0 = Vi- ■ 

C. Transition between two regions 

Now we will characterize the transition between regions (as 
we have seen above, the intersection can be either a simple 
curve in the case of perfect measurements, or a region with 
non empty interior in the uncertain case): 

Lemma 5. For any i £ {1...., n — 1}, consider System 0 
under Control law 0 with Conditions (0© for i,i + lo 
All trajectories with initial conditions in Yj U Vj|i+i enter the 
regular domain Y l+ \. 

Proof: First, we consider i n— 1. We will follow the 
same reasoning as for the previous lemma, applying LaSalle 
theorem on a domain 

D .2 := {£ £ intM.+ \ x < m&x(z(0),Si n )/k, 

s < max(s(0),s in ), y^ < y(£) < t/ 1 } 

for any y^ £ Yj + i. We can show that the functional 
V(£) = y a {D n ) — y(f) is decreasing on Yj whenever u = Dj. 
Similarly, V(£) is also decreasing on Y l+ \ whenever u = 
Dj + 1 . Now under Control law 0, we have shown that V(£) 
is decreasing on the regular domains Yj and Y l+ j . V(£) is 
a regular C 1 function, and can be differentiated along the 
differential inclusion. On the switching domains Yj ij + i, we 
have: 

(dV dV j 

V ® €W {faf& Di) ’fa f & Di+l) f <0 ' 

Thus, V(£) is decreasing on (12- Following the proof of 
Lemma|4]concerning the boundaries, we can deduce that every 
trajectory will reach the boundary y(f) = yf i.e. it will enter 

Y i+1 . 

For i = n — 1, taking any y 1 £ (y n _ 1 , y a (D n )), we can show 
similarly that V(£) is decreasing on O 2 so every trajectory 
will enter the region Y n . ■ 

Following the same proof, we can show that the reverse 
path is not possible, in particular for the last region: 

Lemma 6. Consider System 0 under Control law 0 with 
Conditions 0 for n — 1 ,n. Condition 0 for n — 1, and 
Condition 0- The regular domain Y n is positively invariant. 

D. Global dynamics 

Now we are in a position to present the main result of the 
paper: 

2 or if i = n — 1, Conditions {5) for n — 1, n. Condition {6} for n — 1, and 
Condition 0. 
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Proposition 2. Control law 0 under Conditions 00 
with perfect or uncertain measurements (A1 or A2) globally 
stabilizes System 0 towards the point £, a {D n ). 

Proof: From Lemmas [5] and [6] we can deduce that every 
trajectory will enter the regular domain Y n , and that this 
domain is positively invariant. 

System 0 under a constant control u(t) = D n has two 
non-trivial equilibria (see Proposition0: £ a (A), and £, b {D n ). 
The growth proxy at these two points satisfy y a {D n ) > y n _i 
and yt, (D n ) < y^ (Conditions (|5 1 7 b ) . so there is only one 
equilibrium in Y n : ( If,). Moreover, it is easy to check 

that Y n is in the basin of attraction of £ a (A), therefore all 
trajectories will converge toward this equilibrium. 


V. Implementation of the control law 
A. How to fulfill Conditions 00 

The global stability of the control law is based on Con¬ 
ditions 00. We now wonder how to easily check if these 
conditions hold, or how to choose the dilution rates and/or to 
define the regions in order to fulfill these conditions. In this 
purpose, a graphical approach can be used. 

As an example, we will consider the case where the regions 
are imposed (by technical constraints) and we want to find 
the different dilution rates Di such that Conditions 00 hold. 
Our objective is to globally stabilize the equilibrium point 
£ a (D*) £ Y n , with p(si n ) < D* < pis 0 ) (s° is defined just 
after). 

Let <j>(s) := jp(s)(si n — s), which represents the steady 
state productivity. On [0, Sj n ], (f>(s) admits a maximum for 



Note that we impose D* < pis 0 ) given that for any 
D* > pis 0 ), the same productivity can be achieved with a 
smaller dilution rate, leading to a reduced risk of instability. 

Let us denote s c (y) and Sd(y) the two solutions, for y £ 
(0,</>(s°)), of the equation </>(s) = y, with 0 < s c (y) < s 0 < 
Sd(y) < s in . 

For the Haldane growth rate, we have: 



How to choose Di in one region: For given lower and upper 
bounds y i <y i < cf>{s 0 ), we define: 

D(y.) := min {p(s 0 ), {posdfiyj) and Dfifjf) := iposfiiyfi. 

This can be done analytically or graphically, as shown on 
Figure 0A Whenever we choose A £ {Dffiji), D{y )), we 
have 




Fig. 3. How to choose all the Di (see Section rm Each color represents 
a region (1: blue, 2: green, 3: red, 4: cyan). Top: perfect measurements (Al); 
Bottom: uncertain measurements (A2). 


. s° < S d (y.) < s b {Di), so y b (Di) < y. (Condition 0), 

• sdyf) < s a (A) < so 2/„(A) > Vi (Condition ©). 

so Conditions 0© for i hold. 

If Dfyf) > D(y ) (see Figure 03), it is not possible to 
fulfill the conditions and thus to implement the control law 
with this measurement range. 

How to choose all the Dp The procedure proposed in the 
previous subsection should be repeated for all the regions. We 
can depict two particular cases: 

• For the first region Yj, given y i = 0, we actually impose 

A < p{s in ). 

« For the last region Y n , one should check that Dfiy n _ i) < 
D n = D* < D(y n ). 

This approach is illustrated on Figure 0 

Increasing measurement resolution: We have seen that for 
a given region, it is not always possible to fulfill Conditions 
00. This gives rise to a question: is there any constraint on 
the measurement that guarantees the possibility to implement 
the control law? 

For perfect measurements (Al) with equidistribution, we 
will show that it will always be possible to implement the 
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Fig. 2. How to choose D in one region (see Section [V-AF A: One can choose any D (._ (D(?i). D(y))\ B: it is impossible to fulfill Conditions jnlfil . given 
that D.(y) > O(y)); C: 'H'(y) represents the lower bound for y (limit case). 


control law increasing the number of regions. 

First, we can arbitrarily define the lower bound of the last 
region: y £ (yb(D*),y a (D*)), recalling nonetheless that n 
is unknown. 

Now, we will determine the limit range of a measurement 
region. Let y such that D{y ) = y(si n ). For y . < Vi < y, 
we have < y{si n ) < D{y^ so it is always possible to 

choose a Di in order to fulfill Conditions l(5][6]i. 

For y £ [y, y ], we define (see Figure [2jT): 

i>(y) ■= (</>o s fc o yo s c )(y). 

It gives a lower bound on y : whenever one choose y . > i/j(l ] i ), 
then we have: 

D(y i )>D(m))=m), 

so one can choose any Di £ f(D(y i ),D(y i ) \ and Conditions 
(EEI for i will hold. 

We can now define the mapping A : y i—>■ y — ip(y) on 
[y,y n ], A (y) defines the maximal range of the region with 
upper bound y (in order to be able to find a dilution rate such 
that Conditions (06l hold). A admits a minimum: 

A m := _ min A(y) > 0. 
y6[y>y n l 

y 

Thus, for n > -g 1 - + 1, the regions Yi defined by: 


allow the implementation of the control law. In conclusions, 
whenever the measurement resolution is good enough (i.e. the 
number of regions is high enough), it is always possible to find 


a set of dilution rates Di such that Conditions (EH) hold in 
the case of perfect measurements (Al) with equidistribution. 
For uncertain measurements (A2), the same result holds, but 
the proof is omitted for sake of brevity. 

B. Simulations 

As an example, we consider the anaerobic digestion process, 
where the methane production rate is measured. Parameters, 
given in Table Q] are inspired from ED (considering only 
the methanogenesis step). For uncertain measurements, we 
use discrete time simulation. At each time step tk (with 
At = 0.05d), when £(ffc) is in a switching region Yhj + i, we 
choose randomly the control u(tk) between I), and D l+ i. In 
this case, we perform various simulations for a same initial 
condition. 

TABLE I 

Parameter values used for simulation. 


Parameter 

Value 

p- 

0.74 d -1 

ks 

0.59 g.L —1 

ki 

16.4 g.L— 1 

k 

30 

a 

11 L CH 4 .g —1 

Sin 

30 g.L- 1 


Our objective is to stabilize the equilibrium £ a (D*), with 
D* = 0.47 d 1 (which corresponds to a productivity of 92% 
of the maximal productivity). We first consider the following 
perfect measurement set (with equidistant region): 

yi = h+i = * = 1 »-> n - 1 > 
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(a) Conditions J5171 are fulfilled 
(see Fig. [4) 


(b) Condition 0 -2 does not hold 
(see Fig. [6) 


(c) Condition © -2 does not hold 
(see Fig. [7} 


Fig. 5. Transitions between regions. Ellipses: regular domains; rectangles: switching domains. White regions are transient, while grey regions have a stable 
equilibrium. Ellipses with a thick line are positively invariant. 


with y n = 4 L CH 4 .L 1 .d 1 . With four regions (n = 4), we 
can define a set of dilution rates such that Conditions m 
are fulfilled (see Fig. 0: 

Di = 0.19 d -1 , D 2 = 0.29 d -1 , 

D 3 = 0.4 d^ 1 , 1)4 = 0.47 d _1 . 

For uncertain measurements, we increased each upper 
bound and decreased each lower bound by 10%. It appears 
that the same dilution rates can be chosen. 

Trajectories for various initial conditions are represented in 
the phase portrait for perfect and uncertain measurements, see 
Fig. [4] In accordance with Proposition [2] all the trajectories 
converge towards the set-point. Thus, the transition graph is 
deterministic (there is only one transition from a region to the 
upper one), see Fig. [5a] 

If the number of regions is reduced (three regions only), it is 
not possible in this example to choose dilution rates such that 
Conditions (07]) hold for all i. In this case, some trajectories 
do not converge towards the set-point. Some regions have 
transitions towards the upper region, but also towards the 
lower one. There are sliding modes. This aspect will be further 
discussed in the next subsection. 

C. When conditions are not verified: risk of failure 

We here detail what happens if Conditions 03 are not 
fulfilled, and in particular if there is a risk of washout. This 
point is illustrated by Fig. [6] and Fig. [7] 

First, given the previous analysis of the system, one can 
easily see that only the condition yb{D 1 ) < y , i.e. Di < 
y(si n ) is necessary to prevent a washout, so D\ can be chosen 
with a safety margin in order to avoid such situation. Now, if 
Condition © is not fulfilled for some i > 1, the unstable 


equilibrium fDf) will be located in the region Y,. Thus, 
the region Y, have transitions towards the lower region, and a 
trajectory can stay in the switching domain Tj_iu. Given that 
z = s + kx converges towards Sj n (cf. proof of Lemma ©. 
such trajectory will converge towards the intersection between 
the switching domain and the invariant manifold z = Si n (see 
Fig. [5b] and Fig. ©: 

• For perfect measurements (Al), this gives rise to a sliding 
mode and the convergence towards a singular equilibrium 
point. 

• For uncertain measurements (A2), all the trajectories 
converge towards a line segment. In our simulation, they 
actually also converge towards a singular equilibrium 
point. 

This situation can be detected by the incessant switches 
between two regions. In such case, the dilution rate Di should 
be slightly decreased. 

On the other hand, if Condition © does not hold for some i, 
the stable equilibrium t; a (Di) will be located in the region Yj, 
so some trajectories can converge towards this point instead 
of going to the next region, see Fig. [5c] and Fig. [7] To detect 
such situation, a mean escape time for each region can be 
estimated using model simulations. A trajectory which stay 
much more than the escape time in one region should have 
reached an undesirable equilibrium. In this case, the dilution 
rate Di should be slightly increased. 

In all the cases, the trajectories converge towards a point 
or a line segment. Although it is not desired, this behavior 
is particularly safe (given that there is theoretically no risk 
of washout). Moreover, as explained above, these situations 
can easily be detected and the dilution rates can be changed 
accordingly (manually or through a supervision algorithm). 










1.1 






Fig. 4. Trajectories (magenta lines) with Con trol law 0 for various initial 
conditions in the phase portrait. Conditions 107} are fulfilled, so all the 
trajectories converge towards the set-point (dark circle), see Proposition [2] 
Open circle: washout. Top: perfect measurements (Al); Bottom: uncertain 
measurements (A2). The frontiers are represented by the green dashed lines, 
switching regions are colored in gray. In these gray regions, the system is not 
deterministic. 

VI. Conclusion 

Given the quantized measurements, we were able to design 
(under some conditions) a control based on regions and tran¬ 
sition between regions. These tools are similar to the ones of 
piecewise linear systems, and it is possible to draw a transition 
graph showing all the possible transitions. Moreover, we have 
seen that for some undesirable cases, singular behaviors (slid¬ 
ing modes) are possible on the boundaries between regions. 
We think that this kind of control on domains, and the design 
of the resulting transition graph, is a promising approach, that 
we want to deepen in future works. This approach could be 
generalized to other classical systems, e.g. in mathematical 
ecology. 




Fig. 6. Trajectories with Control law 0 when Conditions {6]>-2 are not 
fulfilled . Same legend as Fig.0] Some trajectories converge towards a singular 
equilibrium point (black diamond) with a sliding mode (see also the associated 
transition graph in Fig. HE}. 
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VII. Appendix 


Proof of Lemma \2\ 

Let us define <p J D (s) := go{s) — h J D (s) for j = a,b. We 

j ^ _ D (sj{D) - s) fi{s) - D 
k/i(s) k/j,'(s) 


<Pd 0) = 


< 4 ( s ) = 

n(s) — D / 1 


li{s) //(s ) 2 


k/i(s) 2 ' 


First, we consider on (0, s). Given that y,(s) is 

increasing and concave on this interval, we get < 0 

on (0 ,s a {D)), and > 0 on (s 0 (£>),s). Moreover, we 

have ipf)(s a (D )) = 0 , so > 0 on ( 0 ,s), which proves 

(i). 


Now we want to determine the sign of p b D (s) on (s, +oo). 
For this purpose, we consider the equation <p b D (s) = 0. By 
replacing /z(s) and its derivative by their analytic expressions, 
this equation becomes: 

S AP1 S 2 - 2 k s s + k s s b {D) = 0. 
ki 

Given that s a (D)sb(D) = kski , the equation <p b jj(s) = 0 has 
only one root s = Sb(D). Moreover, we have: 

lim (fi b ]j(s) = —oo and lim p b £,{s) = —oo. 

s—H-oo 

Given that ^d( s ) i s continuous on (s,+oo), we finally 
conclude that on this interval, <p b D (s) < 0 , i.e. <?d(s) < 
h b D {s). ' □ 


Proof of Lemma [j] 

First, given that <?_d(s) represent the nullcline y(£) = 0, we 
can check that we have y(£) > 0 on (see Figure Q}: 

{£ e | s < s, x < go{s)} 

U {£ (E | s > s, x > 3d(s)}- 

Recalling that h a D (s) and h b D (s ) are respectively the isolines 
y (0 = Va{D) and y(£) = y b (D), Lemma [2] allows to 
conclude that for £ £ such that yb{D) < y(£) < y a {D), 
we have y(Jf) >0. □ 












